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In a recent funded Key Project of National Science Foundation of China, we have planned to do large scale 
simulations of lattice Quantum Chromodynamics, using the parallel supercomputing facilities in mainland China. 
Here I briefly review our plan and progress in recent years. 



1. INTRODUCTION 

Lattice gauge theory (LGT) has been accepted 
as the most powerful tool for investigating non- 
perturbative strong interactions. Chinese physi- 
cists at Institute of High Energy Physics (IHEP), 
Institute of Theoretical Physics (ITP) and Peking 
Univ. (PKU) in Beijing, Nankai Univ. in Tian- 
jin, Chengdu Univ. of Science and Technology 
and Sichuan Univ. in Chengdu, Zhejiang Univ. 
(ZJU) in Hangzhou, and Zhongshan Univ. (ZSU) 
in Guangzhou have been involved in the study of 
LGT since early 80's. However, earlier investiga- 
tions were mainly analytical ones (for a review, 
see Refs.niB ; )■ due to limited computational fa- 
cilities. 

Thanks to the rapid development of high per- 
formance supercomputers in China in late 90's, 
success of the Symanzik improvement program, 
and increasing support from NSF (National Sci- 
ence Foundation), more and more Chinese lattice 
physicists are doing numerical simulations. 

2. COMPUTING FACILITIES 

There have been many commercial supercom- 
puters installed at research institutions and uni- 
versities in China. The Beowulf clusters, however, 
are becoming more and more popular. 

The ZSU group built a cluster|SI of 10 PC type 
computers. The major difference in our comput- 
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ers from one likely to be found in a home or 
business is that each is equipped with two CPUs. 
This allows us to roughly double our processing 
power without the extra overhead cost for extra 
cases, power supplies, network cards, etc. Specif- 
ically, we have installed two Pentium 3 processors 
in each motherboard. 

The Legend Group built a cluster, named 
DeepComp 1800/LSSC-II, which is installed at 
Academy of Mathematics and System Science, 
Beijing, and also open to the Chinese lattice 
community. It consists of Myrinet and 512CPU 
(Xeon-2GHz). Its Rmax and Rpeak are 1.046 
Tflops and 2.048 Tflops respectively. It was the 
number 43 of top 500 supercomputers^] and the 
3rd fastest in Asia in 2002. 



3. NATIONAL SCIENCE FOUNDATION 

The Chinese NSF Committee |S] was estab- 
lished in 1986, and it is an organization directly 
afhliated to the State Council for the manage- 
ment of the National Natural Science Fund. As 
shown in Tab. ^ these years have seen increas- 
ing funding. Each General Project is funded for 
100 - 200K yuan/3 years (1 Yuan=l/8USD); 
Each Project for Distinguished Young Scientists 
is funded for 0.8M yuan/4 years; Each Key 
Project is funded for 1 ~ 2M yuan/4 years. 

Chinese lattice physicists have got a lot of sup- 
port from the NSF, including the fund for Distin- 
guished Young Scientists and the Key project|n|. 
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Table 1 

Approved funding by the Chinese National Science Fundation (in MiUion Yuan). 



Year 


1996 


1997 


1998 


1999 


2000 


2001 


Total 
% of increase /year 


645.8 
26.4 


777.2 
20.4 


888.6 
14.3 


1083.4 
21.9 


1284.3 
18.5 


1568.4 
24.5 



4. PLAN 

In the recently funded Key NSF Project]^, we 
have planned to do large scale simulations of lat- 
tice QCD, using the parallel supercomputing fa- 
cilities in China. We plan to develop new numer- 
ical methods and study the following hot topics: 
new hadrons such as glueballs and hybrid mesons, 
scattering of hadrons, topology of QCD vacuum, 
transition from the quark confinement phase to 
quark-gluon plasma phase, quantum instantons 
and quantum chaos. These investigations will 
provide useful information for the frontier exper- 
iments in particle physics. 

5. RECENT PROGRESS 

5.1. Improved actions and Hamiltonians 

The idea of Symanzik improvement [3| is to re- 
duce the lattice spacing errors by adding new 
terms to the Wilson actions. The pursuit of 
the Symanzik program, combined with tadpole 
improvement has recently led to significant 
progress in LOT, opening the possibility of ap- 
proaching continuum physics on coarse and small 
lattices. Chinese lattice physicists also con- 
tributed actively to the improvement program: 
the development of the next to nearest neighbor 
improved quark actionj^ and improved Hamilto- 
nians for the Quarks [TUI and e:luons|TT|. 

5.2. Glueballs and hybrid mesons 

Lattice computation of the glueball masses 
has a long history. Early investigations 
used the unimproved Wilson gluon act ion [T^ 
[T^ and Kogut-Susskind gluon Hamiltonian jl4ll51 
I16 |. Recent calculations used the tadpole im- 
proved action [r7. 18 19 2 0121) on anisotropic lat- 
tices. In R.ef. |19l2T)] . the IHEP group performed 
a detailed study of the continuum behavior of the 
glueball J = 0, 2, 4 operators and computed the 
corresponding spectrmn. In Ref. 21 , the ZSU 



group performed the calculation of the 0++, 2++ 
and T^ glueball masses, at (3 = 2.2, 2.4, 2.6, 
3.0, and 3.2, using the improved gluon action on 
the anisotropic lattices J7]. In comparison with 
Ref. we have data at weaker coupling and 
larger lattices, so that the extrapolation to the 
continuum limit might be more reliable. Our 
results for the glueball masses are Mg{0~^~^) = 
1737(±52)(±107), Mg(2++) = 2409(±35)(±103) 
and Mg(1+") = 2946(±63)(±142). 

There have been several quenched lattice 
calculations [221 of hybrid meson masses using re- 
spectively (1) the Wilson quark action, (2) the 
SW improved quark action, (3) the NRQCD ac- 
tion, (4) the LBO action, and (5) the improved 
KS quark action. In those simulations, config- 
urations were generated using the Wilson gluon 
action. 

In Ref. |23j . the ZSU group employed both im- 
proved gluon and quark actions on the anisotropic 
lattice, which should have smaller systematic er- 
rors, and should be more efficient in reducing the 
lattice spacing and finite volume effects. We ob- 
tained 2013 ± 26 ± 71 MeV for the mass of the 
1 ^ hybrid meson in the light quark sector, and 
4369 ± 37 ± 99MeV in the charm quark sector; 
the mass splitting between the 1 ^ hybrid me- 
son in the charm quark sector and the spin av- 
eraged S-wave charmonium was estimated to be 
1302 ± 37 ± 99 MeV. 

5.3. QCD at finite chemical potential 

The Hamiltonian formulation of lattice QCD 
doesn't encounter the complex action problem. In 
Refs. |24l25] . the ZSU group developed a Hamilto- 
nian approach at finite chemical potential fi and 
at zero temperature obtained reasonable results 
in the strong coupling regime. A first order chiral 
phase transition was observed at some fic ■ Quan- 
tum chaos and quantum instantons 26 were also 
investigated; These phenomena might be relevant 
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for the rich phase structure. 

5.4. Hadron scattering 

The lattice collaboration in Beijing studied 
scattering of hadrons using tadpole improved 
clover action on coarse anisotropic lattices |^7I28| . 
The scattering lengths of tttt in the 1=2 chan- 
nel and KN in the 1=1 channel were calcu- 
lated within quenched approximation. Compari- 
son with the new experiment and Chiral Pertur- 
bation Theory was made and good agreements 
were found. 

5.5. Chiral properties of overlap quarks 

The ZJU group studied the chiral symmetry 
relation and scaling of the overlap fcrmions|29| 
I3()| . Divergence of quenched chiral logarithm was 
found in very light quark regime. 

5.6. Algorithms 

In Ref.;31,. the ZSU group developed the 
molecular dynamics equations for doing HMC 
simulations with dynamical clover quarks and 
proposed an even-odd preconditioning method for 
the quark determinant. In Ref.[H2]) the algorithm 
was implemented. 

The PKU group developed an algorithm for 
Ginsparg- Wilson quarks[33]. The ZSU group ex- 
tended an algorithm for computing the thermody- 
namical quantities of LGT with Kogut-Susskind 
quarks, including the chiral limit, to QCD|84|. 

The ZSU group developed an algorithm, named 
MC IIamiltonian_35 to study the excited states 
of LGT and tested in the scalar models [S^. 

Some works mentioned above were done with 
the Laval Univ. and Kentucky Univ. groups. 
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